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Abstract 


This  report  summarizes  results  of  stability  analyses  performed  for  a  slender  sharp  cone  with  and 
without  injection  of  air  through  a  porous  strip.  The  free-stream  parameters  are  related  to  the 
experiments  in  the  Caltech’s  T5  shock  tunnel.  The  analysis  is  focused  on  pure  aerodynamic 
effects  in  the  framework  of  perfect  gas  model.  It  was  shown  that  the  injection  leads  to 
destabilization  of  acoustic  modes  in  the  near- fie  Id  relaxation  region.  To  reduce  this  detrimental 
effect  it  was  suggested  to  alter  the  injector  surface  shape  or  use  suction-blowing  of  zero  net 
injection.  However,  stability  computations  showed  that  these  modifications  did  not  improve  the 
injector  performance.  The  porous  wall  effect  on  the  acoustic  instability  was  also  examined.  It 
was  found  that  porosity  stabilizes  the  boundary  layer  flow  with  and  without  injection.  This 
indicates  that  the  injector  performance  can  be  improved,  if  the  porous  layer,  which  is  used  for  the 
injection,  is  protruded  downstream  to  cover  the  near-filed  relaxation  region. 
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Superscripts 


acoustic  admittance  of  a  porous  layer 
local  speed  of  sound 
phase  velocity 
frequency 

cone  length  along  the  symmetry  axis 

Mach  number 

injection  mass  rate 

amplification  factor 

mean- flow  pressure 

pressure 

mass  flow  rate  of  suction/blowing 
eigenfunction 
Reynolds  number 
temperature 

temperature  disturbance 
time 

mean-flow  velocity  along  the  cone  surface 
velocity  components  of  the  disturbance 
coordinates 

streamwise  component  of  the  wavenumber  vector 

circumferential  component  of  the  wavenumber  vector 

boundary-layer  thickness 

displacement  thickness 

specific  heat  ratio 

density 

spatial  growth  rate 
half-angle  of  the  cone 
porosity 

circular  frequency 

sonic  line 
imaginary  part 
upper  boundary-layer  edge 
real  part 
transition  onset 
on  the  wall  surface 
free-stream  quantity 
neutral  point 

dimensional  quantity 
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1.  Introduction 


It  is  well  known  that  laminar-turbulent  transition  on  slender  bodies  in  a  hypersonic  flow  at  small 
angles  of  attack  is  caused  by  amplification  of  acoustic  waves  trapped  in  the  boundary  layer  [1,2]. 
Non-equilibrium  effects  such  as  molecular  vibrations  and  dissociation  can  damp  acoustic 
disturbances  [3,4],  Carbon  dioxide  has  been  found  to  be  well  suited  to  absorb  acoustic  energy  in 
an  enthalpy  range  relevant  to  realistic  applications  [5],  These  findings  point  to  the  potential 
beneficial  application  of  adding  carbon  dioxide  into  boundary-layer  flows  in  order  to  delay 
transition  onset.  Experimentally,  this  is  attempted  by  injecting  CO2  into  the  boundary  layer 
through  the  wall. 

In  Ref.  [6],  three  injection  schemes  were  studied  both  as  passive  trips  with  no  injection  and  as 
active  trips  with  CO2  injection.  The  first  injector  had  four  rows  of  36  orifices.  The  second 
injector  was  derived  from  the  former  one  by  only  keeping  the  fourth  or  most  downstream  row  of 
orifices.  The  third  injector  consisted  of  a  micro-porous  section.  It  was  found  that  the  four-row 
injector  tripped  the  boundary  layer  to  fully  turbulent  values  near  the  last  row  of  orifices.  The 
transition  Reynolds  number  with  the  one-row  injector  was  decreased  by  more  than  50%  as 
compared  to  the  smooth  cone.  Therefore,  while  not  as  efficient  as  the  four-row  injector,  one  row 
of  holes  still  caused  early  transition.  A  porous  injector  was  also  tested  and  did  not  result  in  early 
transition  when  tested  without  injection. 

Due  to  the  tendency  of  the  discrete  jets  to  cause  transition,  further  studies  were  focused  on  the 
porous  injector.  The  numerical  simulations  [7]  predicted  that  transition  should  occur  immediately 
following  the  injection  of  cold  carbon  dioxide.  The  test  cases  with  air  and  nitrogen  as  the  test  gas 
suggested  the  momentum  of  the  injection  plays  a  dominant  role  in  the  amount  of  amplification 
seen  immediately  downstream  due  to  the  inefficient  heating  of  the  injected  CO2.  To  remedy  this, 
the  injected  carbon  dioxide  was  pre-heated,  resulting  in  a  reduction  of  amplification  in  the  post¬ 
injection  region.  Despite  this  reduction,  the  stability  analyses  still  predicted  transition  to  occur 
earlier  as  compared  to  the  case  without  injection.  To  isolate  the  effects  of  non-equilibrium 
processes  on  the  disturbances,  a  second  configuration  involving  a  similar  cone  with  a  longer 
transpiration  interval  was  investigated.  For  this  case,  the  stability  analyses  predicted  a  window  of 
carbon  dioxide  injection  resulting  in  a  reduction  of  amplification  and  thus  a  delay  in  transition.  In 
Ref.  [8],  transition  delays  were  documented  in  shots  with  CO2  injection.  The  data  showed  a 
general  trend  of  increasing  delay  with  injection  rate,  before  a  sharp  drop-off  at  the  highest 
injection  rate.  However,  the  transition  delay  was  essentially  smaller  than  that  observed  in  the 
free-stream  with  CO2. 

The  foregoing  studies  indicate  that  the  injector  performance  essentially  depends  on  the  balance 
between  the  stabilizing  effect  of  carbon  dioxide  and  the  destabilizing  effect  of  injection.  The 
latter  could  be  diminished  by  seeking  for  optimal  distributions  of  the  blow  rate  and/or  shaping  of 
the  injection  region.  To  reach  this  goal  on  economical  basis,  the  following  strategy  was  chosen: 
1)  Perform  broad-band  feasibility  studies  of  different  blowing  configurations  using  the  perfect 
gas  model.  2)  Select  the  most  promising  configurations  and  check  their  robustness  by  numerical 
modeling  of  CO2  injection  including  non-equilibrium  effects.  3)  Design  the  CO2  blowing  system 
and  test  it  in  high-enthalpy  shock  tunnel. 

This  report  is  related  to  Step  1,  where  air  is  injected  to  the  air  flow  and  air  is  treated  as  a  perfect 
gas.  The  analysis  is  focused  on  pure  aerodynamic  effects,  i.e.  absorption  of  acoustic  energy  by 
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C02  is  not  considered.  It  should  be  noted  that  the  perfect-gas  modeling  is  suitable  for  the  near- 
filed  region  where  the  injected  gas  is  relatively  cold.  Further  downstream,  where  the  injected  gas 
is  heated  and  it  can  effectively  damp  acoustic  instabilities,  the  perfect  gas  model  overestimates 
the  disturbance  growth  rates  and,  thereby,  under-predicts  the  injector  efficiency. 

The  following  cases  are  considered:  A)  Baseline  configuration  -  a  5°  half-angle  sharp  cone  at 
zero  angle  of  attack  with  the  normal-wall  blowing  through  a  porous  wall  strip.  This  configuration 
is  similar  to  that  tested  at  Caltech’s  T5  shock  tunnel  [6,8],  B)  Baseline  configuration  with 
shaping  of  the  injection  region;  C)  Baseline  configuration  with  the  normal- wall  suction  and 
blowing  of  zero  net  injection.  The  free-stream  parameters  correspond  to  the  run  2540  in  the  T5 
shock  tunnel,  which  was  documented  in  Ref.  [7].  For  each  case  the  laminar  basic  flow  is 
computed  using  the  in-house  Navier-Stokes  code  HSFlow,  and  stability  analyses  are  performed 
in  the  framework  of  linear  stability  theory  (LST)  using  the  in-house  LST  code.  These  tools  were 
successfully  exploited  for  stability  and  receptivity  studies  of  high-speed  boundary-layer  flows  [9- 
12], 

Fedorov  et  al.  [13]  showed  theoretically  that  a  passive  ultrasonically  absorptive  coating  (UAC) 
may  effectively  stabilize  the  Mack  second  mode.  Typical  UAC  constitutes  a  thin  micro-porous 
layer,  with  average  spacing  between  pores  being  much  smaller  than  the  boundary-layer 
thickness.  The  transition  experiments  of  Rasheed  et  al.  [14],  which  were  performed  on  a  sharp 
cone  in  the  T5  shock  tunnel,  qualitatively  agree  with  the  theoretical  predictions  [13],  Further 
experimental,  theoretical  and  numerical  studies  (see  Rev.  [2]  and  Ref.  [15])  have  demonstrated 
robustness  of  the  UAC  stabilization  concept.  This  motivated  us  to  analyze  the  UAC  effect  on  the 
boundary  layer  instability  in  combination  with  the  gas  blowing.  The  analysis  is  performed  using 
the  LST  computations  with  the  impedance  boundary  conditions  for  disturbances  on  the  porous 
wall.  The  UAC  parameters  are  relevant  to  the  porous  injector  tested  in  the  T5  shock  tunnel  [6]. 
The  UAC  stabilization  effect  is  evaluated  for  the  baseline  configuration  A  with  and  without 
injection. 

2.  Baseline  configuration  and  numerical  approach 

As  a  first  step,  we  consider  the  baseline  configuration  shown  in  Fig.  1.  A  5°  half-angle  sharp 
cone  has  the  region  128  mm  <  x*  <  169  mm,  where  gas  at  the  wall  temperature  is  uniformly 
injected  in  the  normal  wall  direction.  In  our  computations,  the  total  mass  rate  is  ranged  from  zero 
to  13.5  g/s,  where  the  upper  limit  approximately  equals  the  double  mass  rate  of  the  incoming 
boundary-layer  flow. 

The  free-stream  parameters  correspond  to  the  conditions  of  Run  2540  in  the  GALCIT  T5  shock 
tunnel  [7]:  Mach  number =  5.3,  static  temperature  T*.  =  1323.77  K,  static  pressure 
=  21993  Pa,  density  p*^  =  0.05788  kg/m3,  and  wall  temperature  T*  =  293  K.  The  cone 
length  measured  along  the  cone  axis  is  L*  =  1  m.  The  asterisks  denote  dimensional  quantities. 
The  Reynolds  number  based  on  the  cone  length  is  Rexi  =  4.57  •  106 . 

Calculations  of  the  mean  flow  are  performed  using  the  in-house  Navier-Stokes  code  HSFlow  for 
perfect  gas  of  Prandtl  number  Pr  =  0.72,  specific  heat  ratio  7  =  1.4  and  Sutherland  viscosity- 
temperature  dependence,  with  the  Sutherland  constant  being  110.4K.  The  flow  variables  are 
made  nondimensional  using  the  free-stream  parameters  and  the  cone  length  L*  =  1  m.  Pressure 
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is  scaled  using  the  doubled  free-stream  dynamic  pressure  fX.xUx_ ,  where  U*x  is  free-stream 
velocity. 

The  computational  domain  is  a  rectangle  (Fig.  2).  Its  bottom  boundary  lies  on  the  cone  surface, 
and  the  upper  boundary  lies  in  the  free  stream  above  the  cone-induced  shock.  The  boundary 
conditions  on  the  cone  surface  are  the  no-slip  condition  and  the  isothermal  condition  T  =  Tw . 
On  the  inflow  and  upper  boundaries,  the  conditions  correspond  to  the  undisturbed  free  stream. 
On  the  outflow  boundary,  the  unknown  dependent  variables  are  extrapolated  using  the  linear 
approximation. 


Figure  1  A  5°  half-angle  sharp  cone  with  the  injector  -  baseline  configuration  A. 


The  problem  is  solved  numerically  using  the  implicit  second-order  finite-volume  method 
described  in  Ref.  [9],  The  two-dimensional  (axisymmetric)  Navier-Stokes  equations  are 
approximated  by  a  shock-capturing  scheme  that  allows  for  modeling  of  flow  non-uniformities 
induced  by  the  injector.  The  advection  terms  are  approximated  by  the  third-order  WENO  scheme 
[16]  to  decrease  the  numerical  dissipation.  The  code  algorithm  as  well  as  its  implementations  and 
validations  are  discussed  in  [9], 

The  grid  convergence  study  showed  that  the  basic  laminar  flow  can  be  computed  with 
sufficiently  high  accuracy  using  597x649  grid.  This  grid  is  clustered  in  the  direction  normal  to 
the  cone  surface  so  that  the  boundary-layer  region  contains  approximately  50%  of  nodes.  To 
resolve  flow  non-uniformities  in  the  vicinity  of  injector,  the  grid  is  also  clustered  in  the 
streamwise  direction  near  the  injector  boundaries  as  shown  in  Fig.  3. 


Figure  2  Computational  domain  and  grid. 
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Figure  3  Close  up  view  of  the  computational  grid  near  the  injector. 


The  mean  flow  profiles  were  used  for  the  local-parallel  stability  analysis.  The  boundary-layer 
disturbance  is  considered  in  the  wave  form  q(y)exp(iax  +  i(3z  —  iut) ,  where  u  is  angular 
frequency  and  a,  (3  are  streamwise  and  azimuthal  components  of  the  wavenumber  vector,  t  is 
time.  The  vector  q  consists  of  the  eigenfunctions  of  velocity  components  u,v,w,  pressure  p  and 

temperature  T' .  For  each  a; -station  the  linear  stability  equations  are  solved  using  a  4th-order 
Runge-Kutta  scheme  and  a  Gram-Schmidt  orthonormalization  procedure.  For  the  temporal 
stability  problem,  u>(a,/3,x)  is  a  complex  eigenvalue  while  a,  (3  are  real  parameters.  For  the 
spatial  stability  problem,  a(cv,  / 3 ,  x )  is  a  complex  eigenvalue  while  (3,  u  are  real  parameters.  The 
eigenvalues  of  discrete  spectrum  are  calculated  with  the  help  of  a  shooting/Newton-Raphson 
procedure.  The  downstream  amplification  factors  N(x,  u (3)  are  computed  using  solutions  of  the 
spatial  stability  problem 

X 

N(x,uj,/3)=  J  <j(u,/3,x)dx,  (1) 

x0(u,f3) 

where  a  —  —at  is  the  spatial  growth  rate,  and  x0  is  the  initial  point,  which  can  be  the  neutral 
point  on  the  lower  neutral  branch  or  a  certain  fixed  point.  The  subscript  ‘f  denotes  the  imaginary 
part  of  a  complex  quantity  and  V  denotes  the  real  part. 

2.1  Results 

The  mean  flow  fields  of  static  pressure  and  radial  velocity  are  shown  in  Figs.  4a, b  for  the  case  of 
injection  with  the  total  mass  rate  m  =  13.5  g/s  (the  corresponding  nondimensional  mass  flow 
rate  is  pwvw  —  0.0181).  The  close  up  views  of  axial  velocity,  static  pressure  and  temperature 
fields  in  the  injector  vicinity  are  shown  in  Fig.  5,  where  the  dashed  lines  indicate  the  injection 
boundaries.  It  is  seen  that  the  normal-wall  blowing  displaces  the  boundary  layer  and  induces 
compression  waves  emanating  from  the  upstream  boundary  of  the  injector.  A  cold  dead- flow 
layer  is  formed  near  the  wall.  Its  thickness  is  several  times  larger  than  the  shear-layer  thickness. 
Owing  to  the  mixing  process  the  dead-flow  layer  gets  thinner  and  it  is  ultimately  swallowed  by 
the  shear  layer  in  the  mid  stations  (x  «  0.6 ).  This  trend  is  clearly  seen  in  Fig.  6,  where  profiles 
of  the  streamwise  velocity  and  temperature  are  shown  for  various  x  -stations.  Figure  7  shows  the 
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near-wall  layer  thickness  8{x)  determined  by  the  condition  U{8)  =  0.99Ue .  After  a  sharp  rise  in 
the  injection  region,  8(x)  decreases  and  slowly  approaches  the  no-blow  distribution.  These  CFD 
data  demonstrate  that  the  injection  strongly  affects  the  near- wall  flow  not  only  near  the  injector 
but  also  in  a  long  post-injection  region  which  is  called  hereafter  as  a  relaxation  region. 

To  get  insight  into  stability  of  the  near- wall  layer  with  injection,  we  start  with  the  temporal 
stability  analysis  for  the  mean  flow  profiles  in  the  station  x  =  0.298  which  lies  in  the  middle  of 
the  relaxation  region.  The  flow  parameters  are  non-dimensionlized  using  their  values  at  the  upper 
edge  of  near-wall  layer,  the  length  scale  is  the  displacement  thickness  8*  and  the  time  scale  is 
8*  /  U* .  In  the  considered  herein  x  -station,  the  local  Mach  number  is  Me  &  5  and  the  Reynolds 
number  is  Re5,  «  8074 . 
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Figure  4  Basic  flow  fields  for  the  static  pressure  (a),  and  radial  velocity  (b),  the  total  mass  rate 
rh  =  13.5  g/s. 

Using  a  global  search  for  eigenvalues  we  found  seven  unstable  modes  whose  increments  ui^a) 
and  phase  speeds  cr(a)  =  cor(a) / a  are  plotted  in  Figs.  8  and  9.  The  computations  were 
performed  for  plane  waves  of  j3  —  0.  Figure  10  shows  the  pressure  eigenfunctions  p(y)  for 
modes  0,  1,  2  and  3  with  the  wavenumbers  a  corresponding  to  maximal  increments  cu,  for  each 
mode.  The  eigenfunction  shapes  indicate  that  we  are  dealing  with  trapped  acoustic  waves 
propagating  in  a  waveguide  formed  between  the  wall  and  the  relative  sonic  line 
y  =  ya  :  U(ya )  =  cr  —a(ya),  where  a  is  local  speed  of  sound  [2],  In  the  layer  0  <  y  <  ya , 

where  the  mean  flow  is  supersonic  relative  to  the  disturbance  propagating  with  the  phase  velocity 
cr,  the  eigenfunction  pn(y)  of  mode  n  oscillates  versus  y  so  that  the  oscillation  phase  is 
approximately  tt  /  4  +  mr ,  n  —  0, 1, ....  In  the  Mack  terminology,  mode  with  n  =  0  corresponds 
to  the  second  mode,  mode  with  n  =  1  -  to  the  third  mode  etc.  Because  of  relatively  thick  dead- 
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flow  layer  the  wave  guide  is  essentially  thicker  and  the  number  of  unstable  acoustic  modes  is 
larger  than  in  the  no-blowing  case. 

It  turned  out  that  the  phase  speeds  of  acoustic  instabilities  are  quite  small  compared  with  the  no¬ 
blowing  case.  As  shown  in  Fig.  1 1,  the  instabilities  are  observed  in  a  wide  range  of  cr  and  their 
eigenvalues  are  close  to  the  branches  of  slow  acoustic  waves  of  the  continuous  spectrum.  The 
eigenfunctions  of  such  slow  modes  oscillates  outside  the  boundary  layer  with  very  weak 
damping  (Fig.  12).  This  abnormal  behavior  may  lead  to  a  significant  increase  of  receptivity  to 
free-stream  acoustic  disturbances. 


Figure  5  Basic  flow  fields  for  the  axial  velocity,  pressure  and  temperature  in  the  vicinity  of 
injector,  the  total  mass  flow  rate  rh  =  13.5  g/s.  The  dashed  lines  show  the  injection  boundaries. 
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U-velocity 


(a) 


(b) 

Figure  6  Profiles  of  the  streamwise  velocity  U(y)  (a)  and  temperature  T(y)  (b)  at  various  x- 
stations,  the  total  mass  flow  rate  rh  =  13.5  g/s,  the  injection  region  is  shown  by  the  red  line. 
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0.0030 


Figure  7  The  boundary-layer  thickness  in  the  no-blow  case  (black  line)  and  in  the  case  of 
injection  with  m  =  13.5  g/s  (red  line). 


unstable  acoustic  modes  at  x=0.298,  Me=5,  Re  =8074 

0 


streamwise  wavenumber  a 

Figure  8  Growth  rates  cot(a)  of  seven  unstable  modes  at  x  —  0.298 ,  (3  —  0 . 
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Figure  9  Phase  speeds  cr(a )  of  seven  unstable  modes  at  x  =  0.298  ,  (3  =  0 . 
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Figure  10  The  pressure  eigenfunctions  pr(y)  and  p{(y)  for  modes  0,  1,  2  and  3  at  (u>,a) 
corresponding  maximal  ui .  The  red  and  black  lines  show  the  mean  flow  profiles  of  streamwise 
velocity  U(y )  and  temperature  T(y) ,  respectively. 


Figure  11  The  phase  speed  trajectories  (red  lines)  of  unstable  modes  in  the  complex  c  plane,  the 
branches  of  continuous  spectrum  are  shown  by  the  black  lines. 
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Figure  12  Eigenfunctions  pr(y)  and  pt(y)  of  mode  0  (Mack  second  mode)  at  the  phase  speed 
cr  =  0.4  .  The  red  and  black  lines  show  the  mean  flow  profiles  of  streamwise  velocity  U(y )  and 
temperature  T(y) ,  respectively. 

Note  that  maximal  temporal  instability  corresponds  to  mode  2  which  is  the  Mack  fourth  mode 
(Fig.  8).  However  the  spatial  stability  analysis  leads  to  the  different  conclusion.  Figure  13  shows 
that  the  spatial  growth  rates  a(u)  are  maximal  for  mode  0  corresponding  to  the  Mack  second 
mode  -  typical  for  hypersonic  boundary-layer  flows.  This  switch  over  is  associated  with 
abnormal  behavior  of  the  group  velocities  of  the  unstable  modes  considered. 

Figure  14  shows  distributions  of  maximal  growth  rates  amax  (x)  =  max  a(uj,  x)  for  the  five 

unstable  modes  in  the  relaxation  region.  In  the  all  x  -stations  considered,  the  most  unstable  is 
mode  0.  Computations  for  oblique  waves  (/3  ^  0)  showed  that  their  growth  rates  are  smaller 
than  in  the  case  of  (3  =  0  -  typical  for  acoustic  instabilities.  These  findings  allow  us  to  restrict 
further  parametric  studies  to  the  Mack  second  mode  of  (3  =  0 ,  which  is  called  hereafter  as  a 
Mack  mode.  It  should  be  noted  that  the  injection  leads  to  significant  destabilization  of  the  Mack 
mode  in  a  long  region  (compare  the  black  and  orange  lines  in  Fig.  14). 

Distributions  of  the  Mack  mode  growth  rates  (upper  panel)  at  various  frequencies  (lower  panel) 
are  shown  in  Fig.  15  for  the  no-blow  case  (blue  lines)  and  the  case  with  injection  of  m  —  13.5 
g/s  (black  lines).  The  corresponding  distributions  of  N-factors  are  shown  in  Fig.  16.  As 
expected,  the  injection  destabilizes  the  Mack  mode  in  the  relaxation  region  that  gives  a  local 
maximum  of  the  N-factor  envelope  at  x  «  0.3  .  Further  downstream  the  maximal  growth  rates 
remain  higher  than  in  the  no-blow  case.  However  the  width  Ax  of  unstable  regions  is  narrowed 
down  that  leads  to  decreasing  of  the  integral  amplification.  As  a  result  the  N-factor  envelope  is 
lower  than  that  of  the  no-blow  case  in  the  range  0.35  <  x  <  0.75 .  In  the  far-field  region 
x  >  0.6  the  mean  flow  is  almost  parallel  and  the  growth  rates  weakly  depend  on  x .  This  causes 
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a  steep  rise  of  the  N-factor  envelope.  Eventually  this  envelope  crosses  that  of  no-blow  case  at 
x  ft  0.75 . 

The  N  factors  for  the  case  of  rh  =  6.75  g/s  are  shown  by  the  red  lines  in  Fig.  16.  As  expected, 
the  near- field  maximum  of  their  envelope  is  significantly  smaller  than  in  the  case  of  rh  =  13.5 
g/s.  However,  the  relaxation  region  becomes  shorter  and  the  almost  parallel  mean  flow  sets  in 
earlier.  The  envelope  of  steep  slope  is  shifted  upstream  and  the  N-factors  become  larger  than  in 
the  no-blowing  and  rh  =  13.5  g/s  cases  for  x  >  0.5  . 


Figure  13  The  growth  rates  a(oS)  and  phase  speeds  cr{u)  of  seven  unstable  modes  at  x  =  0.298  ,  [3  =  0 . 
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X 

Figure  14  Distributions  of  maximal  growth  rates  along  the  relaxation  region  in  the  case  of 
injection  with  m  =  13.5  g/s;  dashed  lines  show  the  injector  boundaries;  the  orange  line  shows 
the  baseline  case  without  injection;  (3  =  0. 
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Figure  15  Growth  rates  and  frequencies  of  the  Mack  mode  for  the  no-blow  case  (blue  lines)  and 
the  case  with  injection  of  m  =  13.5  g/s  (black  lines). 


Figure  16  N-factors  of  the  Mack  mode  for  the  no-blow  case  (blue  lines),  the  case  of  m  —  13.5 
g/s  (black  lines)  and  the  case  of  m  =  6.75  g/s  (red  lines);  dashed  line  -  Ntr  =  9.2 . 
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Thus,  the  considered  injection  leads  to  destabilization  of  the  near- field  relaxation  region, 
stabilization  of  the  mid-field  relaxation  region,  and  destabilization  of  the  far-field  relaxation 
region  where  the  basic  flow  is  almost  parallel.  The  width  and  location  of  these  regions  as  well  as 
the  level  of  stabilization/destabilization  effect  depend  on  the  injected  mass  rate.  To  estimate  the 
transition  onset  we  use  hereafter  the  critical  N-factor  Nlr  —  9.2 ,  which  has  been  derived  in  Ref. 
[7]  from  comparison  of  the  eN  computations  with  the  transition  measurements  on  a  5  half-angle 
sharp  cone  without  injector,  which  was  tested  in  the  GALCIT  T5  shock  tunnel.  For  sufficiently 
large  rh ,  at  which  the  maximum  of  N-factor  envelope  iVmax  in  the  near-filed  relaxation  region 
is  higher  than  the  critical  N  =  Ntr,  the  transition  onset  point  xtr  is  predicted  to  be  close  to  the 
injector.  Decreasing  of  rh  weakly  affects  xtr  until  _/Vmax  becomes  lower  than  Ntr.  At  this 
moment  the  transition  point  jumps  to  the  x  -station,  which  may  be  even  downstream  from  the 
transition  point  in  the  no-blowing  case;  i.e.,  the  premature  tripping  may  be  switched  to  delay  in 
transition.  Further  decreasing  of  rh  to  a  certain  level  can  move  the  transition  point  upstream. 
Note  that  the  data  in  Fig.  16  are  consistent  with  those  of  Ref.  [7]  where  stability  computations 
were  performed  for  CO2  injection  including  real-gas  effects.  This  convinces  us  that  the  perfect 
gas  modeling  can  be  used  for  a  first-cut  evaluation  of  the  injector  performance. 


3.  Shaping  of  the  injector 

The  foregoing  stability  analysis  suggests  that  the  injector  performance  may  be  improved  by 
shaping  of  the  injection  region.  It  was  suspected  that  decreasing  of  the  injector  surface  slope  can 
partially  compensate  the  displacement  effect  induced  by  the  blowing  and,  presumably,  reduce  the 
instability  growth  in  the  relaxation  region.  This  speculation  was  examined  by  parametric  stability 
studies  of  the  following  configurations. 

The  first  configuration  is  schematically  shown  in  Fig.  17.  The  injection  region  has  a  conical 
shape.  The  rear  part  of  the  original  cone  (x  >  169  mm)  is  shifted  along  the  x  -axis  to  A  b,  and  it 
is  connected  with  the  fore  part  through  a  conical  surface  (red  line)  used  for  injection.  The  angle 
of  the  new  injector  surface  decreases  from  9C  =  5°  to  9C  —  A 9  and  its  x  -length  increases  from 
b  =  41  mm  to  b  +  Ab .  For  this  configuration  stability  computations  were  performed  for  the 
injector  shapes  corresponding  to  A 9  =  1°  and  3°.  The  total  mass  flow  rates  were  rh  =0,  6.75 
and  13.5  g/s  for  each  shape. 

The  second  configuration  has  the  injector  of  cylindrical  shape  shown  in  Fig.  18.  In  this  case 
A 9  —5°  and  the  injector  length  is  41  mm.  Stability  computations  were  performed  at  rh  =  0, 
1.35,  4.05,  6.75,  9.45  and  13.5  g/s. 

The  third  configuration  is  a  combination  of  cylindrical  and  conical  parts  shown  in  Fig.  19.  These 
parts  of  the  injector  surface  are  connected  at  the  mid  point  (xt  +  x2)/2.  Stability  computations 
are  performed  at  rh  —  13.5  g/s. 
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Figure  17  Injector  of  conical  shape. 


Figure  19  Injector  of  cylindrical-conical  shape. 
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3.1  Conical  injectors 


Stability  computations  for  the  conical  injector  shapes  of  A 9  =  1°  and  3°  at  zero  blowing  showed 
that  the  shaping  weakly  affect  N  factors  of  the  Mack  mode.  As  an  example,  Figure  20  compares 
the  distributions  N(x)  for  the  baseline  case  (straight  cone)  with  the  case  of  AO  —  3°.  The 
corresponding  distributions  of  the  boundary-layer  thickness  (Fig.  21)  indicate  that  the  shaping 
without  blowing  produces  a  local  effect  on  the  mean  flow.  For  x  >  0.2  and  x  <  0.128,  the 
distributions  are  close  to  each  other. 
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Figure  20  N-factors  of  the  Mack  mode  for  the  straight  cone  (black  lines)  and  for  the  shape  of 
A 6  —  3°  (red  lines)  at  m  —  0  g/s. 


Figure  21  The  boundary-layer  thickness  for  the  straight  cone  (black  line)  and  for  the  shape  of 
A 0  =  3°  (red  line)  at  rh  =  0  g/s. 
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Now  we  consider  the  conical  injectors  of  A 6  =0°,  1°  and  3°  at  the  fixed  mass  flow  rate  m  = 
13.5  g/s.  Figures  22  and  23  compare  the  basic  flow  fields  for  static  pressure  and  axial  velocity, 
respectively.  As  expected,  the  injector  shaping  reduces  the  pressure  perturbations  induced  by  the 
boundary-layer  displacement.  However,  the  thickness  of  dead-flow  layer  and,  as  a  consequence, 
the  length  of  relaxation  region  increase  with  A 9  (see  distributions  of  6(x)  in  Fig.  24). 
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(C) 

Figure  22  Mean  flow  fields  for  the  static  pressure  at  the  total  mass  flow  rate  rh  =  13.5  g/s;  (a)  - 
A 6  =  0°  (cone),  (b)  -  Ad  =  1° ,  (c)  -  A9  =  3° . 
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Figure  25  shows  the  N  factors  computed  for  the  mean  flows  presented  in  Figs.  22  and  23.  It  is 
seen  that  the  injector  shaping  produces  a  weak  destabilization  effect  in  the  near-field  relaxation 
region  0.2  <  x  <  0.4  and  a  noticeable  stabilization  effect  in  the  mid-field  and  far-field 
relaxation  regions.  Because  the  near-filed  maximum  of  the  N-factor  envelope  remains  to  be 
higher  than  the  critical  level  Ntr ,  these  injectors  may  cause  premature  transition  at  x  «  0.25  . 


Figure  24  Distributions  of  the  boundary- layer  thickness  at  A 9  =  0°  (cone),  A 8  =  1°  and 
A 9  =  3° .  The  total  mass  flow  rate  is  rh  =  13.5  g/s. 


22 

Distribution  A:  Approved  for  public  release;  distribution  is  unlimited. 


20 


Figure  25  N-factors  of  the  Mack  mode  for  the  straight  cone  with  zero  blowing  (blue  lines, 
baseline  case),  and  for  injectors  of  A 6  =  0°  (black  lines),  A 6  =  1°  (red  lines)  and  A 6  =  3° 
(green  lines)  at  m  =  13.5  g/s;  dashed  line  -  Ntr  =9.2. 


3.2  Cylindrical  injector 

Consider  the  cylindrical  injector  depicted  in  Fig.  18.  Close  up  views  of  the  pressure  and  axial 
velocity  fields  are  shown  in  Figs.  26  and  27  for  the  case  of  m  =  13.5  g/s.  Comparing  these 
fields  with  the  baseline  case  (Figs.  22  and  23)  we  conclude  that  the  shaping  reduces  the  pressure 
wave  induced  by  the  blowing  -  consistent  with  the  original  expectation.  The  dead-flow  region, 
however,  becomes  essentially  thicker. 
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Figure  26  Static  pressure  near  the  injector  of  A 6  =  5° ,  0.11  <  x  <  0.3 ,  m  =  13.5  g/s. 
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Figure  28  shows  the  wall-pressure  distributions  Pw(x)  of  the  mean  flow  for  different  injection 
rates  rh  from  0  to  13.5  g/s.  The  dashed  lines  correspond  to  the  injection  boundaries.  For  a 
straight  cone  without  injection  the  baseline  distribution  is  approximately  uniform  (the  solid  black 
line).  The  shaping  without  blowing  (the  red  line)  leads  to  the  pressure  drop  on  the  injector 
surface.  This  behavior  is  consistent  with  the  slender  body  theory  for  inviscid  flows.  The  blowing 
causes  a  sharp  increase  of  Pw  near  the  upstream  boundary  of  the  injector.  Then  Pw  decreases  to 
the  level,  which  is  lower  than  in  the  baseline  case.  Further  downstream  (in  the  relaxation  region) 
the  wall  pressure  is  almost  constant.  The  length  of  Pw  plateau  increases  while  its  level  decreases 
as  the  injection  mass  flow  rate  increases. 

Figure  29  shows  the  corresponding  distributions  of  the  boundary-layer  thickness  6(x) .  In  the  no¬ 
blowing  case,  the  shaping  effect  is  local  and  relatively  weak.  As  the  injection  rate  increases,  the 
boundary- layer  thickness  quickly  increases  in  the  injection  region.  Its  maximum  at  rh  =  13.5  g/s 
is  an  order  of  magnitude  larger  than  in  the  baseline  case.  Further  downstream  the  boundary-layer 
thickness  slowly  decreases  approaching  the  unperturbed  (baseline)  distribution. 

These  mean-flow  changes  affect  the  maximal  (vs.  frequency)  growth  rates  of  the  Mack  mode  as 
shown  in  Fig.  30.  The  shaping  itself  (red  line  for  the  case  of  rh  =0)  weakly  influences  on  the 
distribution  of  crmax .  As  rh  increases,  a  sharp  peak  of  <rmax  (x)  is  formed  near  the  upstream 
boundary  of  injection  a:  «  0.13  .  Another  maximum  is  observed  in  the  relaxation  region.  Because 
its  location  is  shifted  downstream  with  rh ,  the  injection  affects  the  N-factor  distributions  in  a 
nontrivial  manner  (Fig.  31).  The  first  maximum  of  the  N-factor  envelope  (which  is  located  in  the 
near-field  relaxation  region)  increases  with  rh ,  while  the  region  of  steep  envelope  (which  is 
located  in  the  far-field  relaxation  region)  is  shifted  downstream.  Because  the  N-factor  envelopes 
are  non-monotonic  vs.  x ,  the  eN  method  predicts  different  transition  onsets  depending  on  the 
choice  of  critical  N.  For  Ntr  —  9.2  (horizontal  dashed  line),  the  most  promising  cases  are  related 
to  the  injection  with  rh  =  6.75  —  9.45  g/s. 

Figure  32  compares  the  N-factors  for  the  cases  of  A 6  —  0°  (cone  without  shaping)  and 
A 9  =  5°  (cylindrical  injector)  at  rh  =  6.75  g/s  and  13.5  g/s,  respectively.  Although  the 
cylindrical  injector  works  better  in  the  mid-field  and  far-field  relaxation  regions,  its  performance 
is  not  improved  in  the  most  critical  near-filed  region  where  the  N  factor  envelope  has  a  local 
maximum. 

Thus,  the  eN  analysis  indicates  that  the  considered  herein  shaping  does  not  prevent  from  early 
transition  in  the  near-field  relaxation  region  at  sufficiently  large  injection  rates  (e.g.,  the  lower 
panel  in  Fig.  32).  For  moderate  rh ,  at  which  N-factors  in  the  near-field  region  are  below  the 
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critical  level  (e.g.,  the  upper  panel  in  Fig.  32),  the  shaping  produces  a  significant  stabilization 
effect  in  the  mid- field  and  far- field  relaxation  regions. 


Figure  28  The  mean-flow  wall-pressure  distributions  in  the  region  0.1  <  x  <  0.5  at  the  total 
mass  rate  from  0  to  13.5  g/s;  dashed  lines  show  the  injection  boundaries. 


Figure  29  Distributions  of  the  boundary-layer  thickness  at  the  total  mass  rate  from  0  to  13.5  g/s; 
the  cylindrical  injector  of  A 8  =  5° . 
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Figure  30  Maximal  (vs.  frequency)  growth  rates  of  the  Mack  mode.  Black  line  -  baseline  case  of 
straight  cone  without  injection.  Colored  lines  -  cylindrical  injector  (A 8  =  5°)  with  different 
mass  flow  rates  from  0  to  m0  =  13.5  g/s.  Dashed  lines  show  the  boundaries  of  injection.  crmax  is 

made  nondimensional  using  the  Blasius  length  scale  yji'lx*  /  U*e  . 


x 

Figure  31  N-factors  of  the  Mack  mode  for  the  cylindrical  injector  (A 6  =  5°)  at  different 
injection  mass  flow  rates  from  0  g/s  to  rh0  —  13.5  g/s;  dashed  line  -  Nlr  —  9.2 . 
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Figure  32  N-factors  of  the  Mack  mode  for  the  straight  cone  with  zero  blowing  (blue  lines),  for 
the  injector  of  A 6  =  0°  (black  lines)  and  for  the  cylindrical  injector  of  A 9  =  5°  (red  lines).  For 
the  upper  plot  the  injected  mass  flow  rate  is  to  =  6.75  g/s,  and  for  the  lower  plot  to  =  13.5  g/s; 
dashed  line  -  Ntr  —  9.2 . 


27 

Distribution  A:  Approved  for  public  release;  distribution  is  unlimited. 


3.3  Cylindrical-conical  injector 


Consider  the  cylindrical-conical  (CC)  injector  of  the  shape  schematically  shown  in  Fig.  19.  The 
gas  blowing  is  performed  in  the  region  from  x1  =  0.128  to  x2  =  0.169  .  It  was  expected  that  this 
configuration  helps  to  reduce  the  thickness  of  dead- flow  layer  downstream  from  the  injector.  It 
was  also  expected  that  acceleration  of  the  near-wall  flow  in  the  vicinity  of  x  =  x2  may  help  to 
stabilize  the  boundary  layer  in  the  near- field  relaxation  region  and,  thereby,  improve  the  injector 
performance.  To  examine  these  presumptions  we  have  computed  the  mean  flow  and  the  second¬ 
mode  growth  rates  for  the  case  of  the  dimensionless  mass  flow  rate  pwvw  =  0.021  that 
approximately  corresponds  to  the  total  mass  flow  rate  rh  =  13.5  g/s.  Close-up  views  of  the 
pressure  and  axial  velocity  fields  are  shown  in  Figs.  33  and  34,  respectively.  Comparing  these 
fields  with  the  case  of  cylindrical  injector  (Figs.  26  and  27)  we  conclude  that  the  CC  shape  leads 
to  the  significant  increase  of  pressure  perturbation,  while  the  boundary-layer  thickness  is  reduced 
in  the  relaxation  region  (see  Fig.  35). 

The  nondimensional  maximal  (vs.  frequency)  growth  rates  of  the  Mack  mode  are  shown  in  Fig. 
36.  Despite  significant  differences  in  the  mean  flow  characteristics,  ermax  (x)  for  the  CC  injector 
is  not  so  far  from  the  case  of  cylindrical  injector  (compare  the  green  and  red  lines).  As  shown  in 
Fig.  37,  the  dimensional  growth  rates  at  various  fixed  frequencies  are  affected  in  the  two 
different  ways.  The  CC  shape  leads  to  the  reduction  of  maximal  <r*(x)  and  the  increase  of  the 
unstable  region  length.  As  a  result,  the  N  factor  envelope  for  the  CC  injector  is  higher  than  in  the 
case  of  cylindrical  injector.  Thus,  we  conclude  that  the  CC  shape  considered  does  not  improve 
the  injector  performance  at  rh  =  13.5  g/s. 
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Figure  33  Static  pressure  in  the  region  0.11  <  x  <  0.3 ,  rh  =  13.5  g/s. 
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Figure  35  Distributions  of  the  boundary-layer  thickness  for  the  cylindrical  injector  and  the 
cylindrical-conical  injector,  the  total  mass  rate  is  rh  =  13.5  g/s. 


Figure  36  Maximal  (vs.  frequency)  growth  rates  of  the  Mack  mode.  Black  line  -  baseline  case  of 
straight  cone  with  zero  injection;  red  line  -  cylindrical  injector,  green  line  -  cylindrical-conical 
injector  at  rh  =  13.5  g/s.  Dashed  lines  show  the  boundaries  of  injection.  crmax  is  made 

nondimensional  using  the  Blasius  length  scale  Jv  y/u;. 
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Figure  37  Growth  rates  and  N-factors  of  the  Mack  mode  for  the  cylindrical  injector  (red  lines) 
and  the  cylindrical-conical  injector  (green  lines)  at  m  =  13.5  g/s. 
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4.  Suction-blowing  of  zero  total  mass  injection 

Another  way  to  improve  the  injector  performance  could  be  a  combination  of  normal-wall  suction 
and  blowing  with  zero  net  injection  as  schematically  shown  in  Fig.  38.  It  is  expected  that  the 
preliminary  suction  of  the  incoming  boundary-layer  flow  could  partially  compensate  negative 
effects  produced  by  the  subsequent  blowing. 


This  expectation  is  examined  by  stability  computations  for  the  following  distribution  of  suction¬ 
blowing.  In  the  region  —l*  /  2  <x*  —  x*0  <1*  /  2,  where  x*Q  =  212.5  mm  is  the  injector  center 

and  f  =  41  mm,  the  dimensionless  mass  flow  rate  q  =  pwvw  is  specified  as 

q(ip)  =  q0  sin<y?,  <p  =  (x*  —  Xg)2n  /l*,  —  7r  <  cp  <  tv  .  (2) 

Computations  were  performed  for  the  total  mass  rate  rh+  =  6.75  g/s  in  the  blowing  region 

Xq  <  x*  <1*  /  2  (a  half  of  the  maximal  injection  mass  rate  considered  in  the  previous  sections). 
The  mean  flow  field  for  static  pressure,  axial  velocity  and  temperature  are  shown  in  Fig.  39. 
Figure  40  compares  the  wall  pressure  distributions  for  the  baseline  (no-blow)  case,  the  uniform 
blowing  case  of  m  —6.75  g/s  and  the  suction-blowing  case  of  rh+  =6.75  g/s.  It  is  seen  that  in 

the  suction-blowing  case  the  relaxation  region  is  shorter  and  the  dead-flow  layer  is  thinner  than 
in  the  case  of  uniform  blowing.  However  the  local  perturbations  in  the  vicinity  of  the  suction¬ 
blowing  system  are  relatively  large. 

Figure  41  shows  the  N  factors  corresponding  to  the  cases  depicted  in  Fig.  40.  As  contrasted  to 
the  original  expectation,  the  suction-blowing  system  destabilizes  the  flow  compared  with  the 
uniform  blowing  case.  Thus,  the  preliminary  suction  does  not  improve  the  injector  performance. 
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Figure  39  Basic  flow  fields  for  the  static  pressure  (a),  axial  velocity  (b),  and  temperature  (c)  at 
rh+  =6.75  g/s. 


Figure  40  The  wall  pressure  distributions  for  the  baseline  (no  blow)  case  (black  line),  the 
uniform  blowing  of  m  =6.75  g/s  (blue  line)  and  the  suction-blowing  of  rh+  =6.75  g/s  (red 

line). 
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Figure  41  N-factors  of  the  Mack  mode  for  the  baseline  (no-blow)  case  (blue  lines),  the  uniform 
blowing  case  of  rh  =6.75  g/s  (black  line)  and  the  suction-blowing  case  of  rh+  =6.75  g/s  (red 

line);  dashed  line  -  Ntr  =  9.2 . 

5.  Porous  wall  effect 

In  this  section  we  consider  the  porous  wall  effect  on  the  Mack  mode  instability.  In  accord  with 
the  model  developed  in  Ref.  [13]  the  spatial  stability  problem  is  solved  with  the  impedance 
boundary  conditions  for  disturbances  on  the  porous  wall 

%  =  A'P„,  ,  (3) 

where  v  is  the  vertical  velocity  disturbance,  p  is  the  pressure  disturbance.  The  acoustic 
admittance  A  depends  on  the  porous  layer  microstructure.  The  porous  wall  of  the  injector 
considered  herein  was  fabricated  from  sintered  316L  stainless  steel  with  the  media  grade  of  10 
pm.  This  wall  has  a  random  porosity  of  the  average  pore  size  d*  &  10  pm.  The  porous  layer 
thickness  h*  =1.6  mm  is  much  larger  than  d* .  The  top  view  of  the  injector  surface  is  shown  in 
Fig.  42.  For  the  first-cut  evaluation  of  the  porous  wall  effect,  we  assume  that  the  acoustic 
admittance  of  the  actual  porous  layer  is  close  to  that  of  a  porous  layer  with  regular 
microstructure.  Namely,  we  consider  a  porous  wall  with  equally  spaced  vertical  cylindrical 
micro-holes.  Its  admittance  is  expressed  in  the  analytical  form  derived  in  Ref.  [13].  For  /i>  d, 
we  get 

A  =  -<t>!Z„,  (4) 
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where  (f>  is  porosity.  The  characteristic  impedance  Z0  is  expressed  in  terms  of  the  complex 
dynamic  density  and  the  complex  compressibility,  which  are  calculated  using  analytical  solutions 
for  acoustic  disturbances  propagating  within  a  long  cylindrical  pore  of  radius  r* . 

In  the  following  stability  computations,  the  pore  radius  is  r*  =  5  pm  (a  half  of  the  media  grade), 
and  porosity  <fi  =  0.2  corresponds  to  the  pore  spacing  Ar* .  The  porous  coating  covers  the  cone 

surface  from  the  injector  to  the  cone  base.  The  free-stream  parameters  correspond  to  Run  2540  - 
they  are  specified  in  Section  2. 


Figure  42  Magnified  image  of  the  porous  surface  for  the  injector  tested  in  the  T5  shock  tunnel. 

5.1  Porous  coating  for  the  baseline  configuration 

First  we  consider  the  baseline  configuration  A  without  injection.  The  growth  rate  distributions 
c t*(x )  of  the  Mack  mode  at  various  frequencies  are  shown  in  Fig.  43.  The  corresponding  N- 
factors  are  plotted  in  Fig.  44.  It  is  seen  that  the  porous  wall  produces  a  strong  stabilization  effect 
for  the  all  frequencies  considered.  In  the  porous  wall  case,  the  N-factor  envelope  lies  well  below 
the  critical  level  Ntr  =9.2  suggesting  that  the  boundary  layer  flow  should  be  laminar  all  of  the 
way  to  the  end  of  the  cone. 

Our  estimates  show  that  the  pore  radius  is  small  compared  with  the  viscous  Stokes  layer  forming 
on  the  pore  surface  for  the  all  frequencies  considered.  Due  to  viscous  dissipation  the  boundary- 
layer  disturbances  quickly  vanish  inside  each  pore  that  provides  effective  absorption  of  the 
acoustic  energy  and,  thereby,  stabilization  of  the  boundary-layer  flow. 

The  foregoing  model  assumes  that: 

.  The  number  of  pores  per  the  instability  wavelength  ( npor )  is  large  (this  allows  us  to 
average  the  boundary  condition  over  the  porous  surface) 
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•  Roughness  of  the  pore  surface  is  negligible  -  the  surface  is  treated  as  aerodynamically 
smooth 

To  check  the  first  assumption  we  have  calculated  npor  (x)  using  the  wavelength  distribution 
X*(x)  for  the  most  unstable  (vs.  frequency)  waves.  Figure  45  shows  that  npor  >  100 

downstream  from  the  injector  (x  >  0.2).  To  check  the  second  assumption  we  have  calculated 
the  roughness  Reynolds  number  [17] 


^ ekk 


U*(k*)k* 

v\k*) 


r 

u\t;Y  k 


K+K'k*, 


(5) 


where  k*  is  the  roughness  height,  which  is  estimated  as  k*  m  d*  =  10  pm.  Figure  46  shows  that 
in  the  region  0.1  <  x  <  1  the  roughness  Reynolds  number  is  essentially  smaller  than  the 
conservative  criterion  Rett  =  25  for  distributed  roughness  [18].  Thus,  the  both  assumptions  are 
valid  for  the  porous  layer  microstructure  considered. 

Now  we  consider  the  case  of  injection  with  rh  =  13.5  g/s.  As  in  the  case  considered  in  Section 
2,  the  gas  blowing  is  performed  in  the  region  128  mm  <  x*  <  169  mm .  Further  downstream 
the  cone  surface  is  assumed  to  be  covered  by  the  passive  porous  layer.  The  stability 
computations  were  conducted  in  the  range  of  0.2  rri  <  x*  <  1  in  -  downstream  from  the 
injector.  Figure  47  compares  the  growth  rate  distributions  cr*(x )  for  the  case  with  the  porous 
coating  (red  lines)  and  without  it  (black  lines).  The  corresponding  N-factors  are  compared  in  Fig. 
48.  As  expected  the  porous  coating  leads  to  significant  stabilization  of  the  Mack  mode. 


X 


Figure  43  Growth  rates  cr*(x)  of  the  Mack  mode  at  various  frequencies;  black  lines  -  solid  wall, 
blue  lines  -  porous  wall;  there  is  no  injection. 
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Figure  44  N-factors  of  the  Mack  mode  at  various  frequencies;  black  lines  -  solid  wall,  blue  lines 
-  porous  wall;  there  is  no  injection. 


Figure  45  The  number  of  pores  per  the  instability  wavelength  npor(x ) . 
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Figure  46  The  Reynolds  number  distribution  Rett  (x)  at  k*  —  10  jam. 
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Figure  47  The  Mack  mode  growth  rates  a*(x )  at  various  frequencies  in  the  case  of  injection 
with  rh  =  13.5  g/s;  black  lines  -  solid  wall,  red  lines  -  wall  with  the  porous  coating. 
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Figure  48  N-factors  of  the  Mack  mode  at  various  frequencies  in  the  case  of  injection  with 
m  =  13.5  g/s;  black  lines  -  solid  wall,  red  lines  -  wall  with  the  porous  coating. 


This  numerical  example  indicates  that  the  CO2  injector  performance  can  be  improved  using  the 
following  approach.  The  porous  layer,  which  is  used  for  the  injection,  should  be  protruded 
downstream  to  the  near-filed  relaxation  region  as  schematically  shown  in  Fig.  49.  The  protruded 
portion  is  placed  on  the  solid  backup  and  it  works  as  a  passive  porous  coating  (there  is  no  suction 
or  blowing  through  this  coating).  In  the  case  considered,  the  protruded  portion  should  cover  the 
cone  surface  from  x*  =  169  mm  to  approximately  x*  «  300  mm.  This  mat  help  to  reduce  the 
first  maximum  of  the  N-factor  envelope  and,  thereby,  avoid  premature  transition.  Further 
downstream,  the  CO2  stabilization  effect  is  expected  to  be  sufficiently  strong  to  keep  the 
boundary-layer  flow  laminar. 
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Figure  49  Schematics  of  the  injector  with  UAC.  The  porous  layer  is  protruded  downstream  from 
the  injection  region  to  cover  the  near- filed  relaxation  region.  The  protruded  part  is  placed  on  a 
solid  backup 
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6.  Conclusions 


Stability  analyses  of  high-speed  boundary-layer  flow  past  a  5°  half  angle  sharp  cone  with  the 
wall-normal  injection  of  gas  through  a  porous  strip  were  performed  using  Navier-Stokes 
solutions  for  the  mean  flow  and  the  linear  stability  theory.  The  cone  model  and  free-stream 
parameters  were  chosen  to  be  similar  to  the  experiments,  which  were  carried  out  at  Caltech’s  T5 
shock  tunnel  to  investigate  the  effect  of  CO2  injection  on  laminar- turbulent  transition. 

It  was  found  that  the  injection  induces  a  cold  dead- flow  layer,  which  is  slowly  swallowed  by  the 
shear  layer  in  the  downstream  relaxation  region.  The  near-wall  flow  behaves  as  a  wave  guide 
which  can  support  several  unstable  modes  of  acoustic  type.  The  spatial  stability  analysis  showed 
that  the  most  unstable  mode  corresponds  to  the  Mack  second  mode.  The  phase  speeds  of  this 
instability  are  close  to  the  phase  speeds  of  slow  acoustic  waves,  and  its  frequencies  are  several 
times  smaller  than  in  the  no  injection  case.  This  may  lead  to  significant  increase  of  receptivity  to 
free-stream  acoustic  noise. 

The  eN  computations  showed  that  the  injection  leads  to  destabilization  of  the  near- field  relaxation 
region,  stabilization  of  the  mid-field  relaxation  region,  and  destabilization  of  the  far-field 
relaxation  region  where  the  basic  flow  is  almost  parallel.  The  width  and  location  of  these  regions 
as  well  as  the  level  of  stabilization/destabilization  effect  essentially  depend  on  the  injected  mass 
flow  rate.  These  results  are  consistent  with  the  stability  computations  of  Wagnild  et  al.  [7],  which 
were  carried  out  for  CO2  injection  including  the  real-gas  effects.  This  convinces  us  that  the 
perfect  gas  modeling  can  be  used  for  a  first-cut  evaluation  of  the  injector  performance. 

In  order  to  improve  the  injector  performance  it  was  suggested  to  decrease  the  injector  surface 
slope.  It  was  suspected  that  this  shaping  could  partially  compensate  the  displacement  effect 
induced  by  the  injection  and  thereby  reduce  the  instability  growth  in  the  near- field  relaxation 
region.  However,  the  eN  parametric  studies  predicted  that  this  shaping  does  not  stabilize  the  near- 
field  flow  and,  therefore,  it  does  not  prevent  from  early  transition  at  sufficiently  large  injection 
rates.  For  relatively  small  rh  at  which  the  N-factors  in  the  near-field  relaxation  region  are  below 
the  critical  level,  the  shaping  produces  a  significant  stabilization  effect  in  the  mid-  and  far-field 
relaxation  regions.  It  is  expected  that  this  stabilization  will  be  enhanced  by  the  injection  of  CO2. 

Stability  computations  were  also  conducted  for  the  injector  of  cylindrical-conical  (CC)  shape  at 
the  total  injection  rate  rh  =  13.5  g/s.  It  was  found  that  the  CC  shape  leads  to  the  reduction  of 
maximal  growth  rates  and  the  increase  of  the  unstable  region  length.  As  a  result,  the  N  factor 
envelope  for  the  CC  injector  is  higher  than  in  the  case  of  cylindrical  injector  -  the  CC  shape 
considered  does  not  improve  the  injector  performance. 

Apparently,  the  foregoing  shaping  is  not  exhaustive.  Additional  parametric  studies  are  needed  to 
identify  an  optimal  shape  at  which  the  injector  does  not  cause  premature  transition  in  the  near- 
field  relaxation  region. 

It  was  suspected,  that  the  injector  performance  could  be  improved  using  a  combination  of 
normal-wall  suction  and  blowing  with  zero  net  injection.  It  was  assumed  that  the  preliminary 
suction  of  the  incoming  boundary-layer  flow  could  partially  compensate  negative  effects 
produced  by  the  subsequent  blowing.  However,  the  eN  computations  showed  that  the  suction¬ 
blowing  system  destabilizes  the  flow  in  the  whole  relaxation  region  compared  with  the  case  of 
blowing  at  rh  =  rh+  =  6.75  g/s. 
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The  porous  wall  effect  on  the  acoustic  instability  was  also  examined  using  the  linear  stability 
theory  with  the  impedance  boundary  conditions  on  the  porous  surface.  It  was  found  that  porosity 
strongly  stabilizes  the  boundary  layer  flow  with  and  without  injection.  This  suggests  that  the 
injector  performance  can  be  improved,  if  the  porous  layer,  which  is  used  for  the  injection,  is 
protruded  downstream  to  cover  the  near-filed  relaxation  region.  This  will  help  to  reduce  the  first 
maximum  of  the  N-factor  envelope  and,  thereby,  avoid  premature  transition.  Further 
downstream,  the  CO2  stabilization  effect  is  expected  to  be  sufficiently  strong  to  keep  the 
boundary-layer  flow  laminar. 

It  should  be  noted  that  the  local  injection  induces  spatial  non-uniformities  of  the  basic  flow, 
which  can  alter  the  instability  growth  rates  via  nonparallel  effects.  Because  our  stability  analysis 
does  not  account  for  these  effects,  the  foregoing  results  should  be  taken  with  great  caution.  The 
PSE  analysis  and  direct  numerical  simulations  are  needed  to  verify  the  eN  predictions  reported 
herein.  Note  also  that  the  perfect  gas  model  does  not  account  for  the  CO2  stabilization  effect 
owing  to  absorption  of  acoustic  energy.  This  model  is  suitable  for  the  near-filed  region  where  the 
injected  gas  is  relatively  cold  and  the  absorption  effect  is  weak.  Further  downstream,  where  the 
injected  flow  is  heated,  the  perfect  gas  model  overestimates  the  disturbance  growth  rates  and, 
presumably,  under-predicts  the  injector  performance.  Nevertheless,  the  perfect  gas  modeling 
provides  a  good  launching  pad  for  further  stability  analysis  including  the  real  gas  effects. 
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